function res = wishart(V, W, v)

D = size(W,1);
denom = 2^(v*D/2) * gammap(v, D) * (det(W))^(v/2);

res = (1/denom) * (det(V))^((v-D-1)/2) * exp(-1/2*trace(inv(W)*V));
